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Abstract: We discuss some properties of recently proposed models of a ghost-free gravity. 
For this purpose we study solutions of linearized gravitational equations in the framework of 
such a theory. We mainly focus on the version of the ghost-free theory with the exponential 
modification exp(D/^ 2 )D _1 of the free propagator. The following three problems are dis¬ 
cussed: (i) Gravitational field of a point mass; (ii) Penrose limit of a point source boosted 
to the speed of light; and (iii) Spherical gravitational collapse of null fluid. For the first 
problem we demonstrate that it can be solved by using the method of heat kernels and obtain 
a solution in a spacetime with arbitrary number of dimensions. For the second problem we 
also find the corresponding gyraton-type solutions of the ghost-free gravitational equations 
for any number of dimensions. For the third problem we obtain solutions for the gravitational 
field for the collapse of both ’’thin” and ’’thick” spherical null shells. We demonstrate how the 
ghost-free modification of the gravitational equations regularize the solutions of the linearized 
Einstein equations and smooth out their singularities. 
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1. Introduction 

It is widely believed that the quantum gravity will ’’cure” a decease of the classical general 
relativity, its singularities. In particular, in the domain of a spacetime near singularities, where 
the curvature becomes large, the Einstein equations should be modified. Such a modification, 
which, for example, is dictated by the string theory, should include additional terms in the 
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gravitational effective action, that are both, higher in curvature and in its derivatives. It was 
proposed many different modifications of the Einstein theory of the general relativity, that, 
in particular, change its infrared and ultraviolet behavior (see, e.g. a review [1]). In this 
paper we discuss a special class of the modified gravity theory with higher derivatives. It is 
instructive to check at first the effects of such a modification in the linearized version of the 
corresponding theory. It is well known that addition of quadratic in the curvature corrections 
modifies the standard Laplace equation for the gravitational potential <p in the Newtonian 
approximation, which takes, for example, the following form 

{l 2 A +1) A (f = 47T/9 . (1.1) 

It is easy to show that such a modified equation for a point mass has a decreasing at the 
infinity solution, which remains finite at r = 0. The parameter l, which is determined by the 
coupling constant for the quadratic in the curvature term in the effective action, provides a 
UV cut-off in the regularized solution (see e.g. [2, 3]). A general analysis of the Newtonian 
singularities in higher derivative gravity models was recently performed in [4], 

A connected problem is a possibility of mini-black hole production in the gravitational 
collapse of a small mass. To study such a problem one may consider it first in the linearized 
version of the theory. If a corresponding solution is regular and its curvature for a small 
mass M is uniformly small in the whole spacetime, one might conclude that in this regime 
higher in curvature corrections are small and can be neglected. In such a case one can expect 
that for the small enough mass M a black hole is not formed. In other words, in such a 
theory there exists a mass gap for the back hole formation. Long time ago in the paper [3] it 
was demonstrated that the gravitational theory with quadratic in the curvature terms in the 
action possesses this property. 

However, in a general case, a propagator in a theory with higher derivatives contains 
extra poles, reflecting the existence of the additional to the gravitons degrees of freedom. As 
a result, the corresponding theory with higher derivatives usually contains ghosts [2, 5, 6]. 
In order to avoid this problem a new type of the modification of the Einstein theory, called 
ghost-free gravity, was proposed [7, 8, 9, 10, 11]. A similar model was proposed earlier in [12] 
(for general discussion see [13, 14]). 

The main idea of this approach is to consider a theory, which is non-local in derivatives. 
Suppose that in the linearized version of such a theory the box operator □ is modified and 
takes the form a(D)D, where a(z) is an entire function of the complex variable z. Then the 
propagator for such a theory does not have additional poles, different from the original pole, 
describing propagation of the gravitons. There exists a variety of entire functions. It is suffi¬ 
cient to use the function a(z) to be an exponent of the polynomial in z. The simplest choice is 
when this polynomial is a linear function, so that the modified box-operator is exp(—□//r 2 )D, 
where //, which has the dimensionality of the mass, is the UV cut-off parameter. One may 
hope that such cut-off regularizes singularities inside black holes [15, 16] and in the big-bang 
cosmology [17, 18, 19, 20, 16]. It should be emphasized that non-local operators of a similar 
form and their properties were considered long time ago in the papers [21, 22, 23, 24, 25, 26]. 
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In the present paper we study solutions of the ghost-free gravitational equations in the 
linearized approximation. In the section 2 we demonstrate that for any number of space- 
time dimensions D a solution for the Newtonian potential for a point mass can be obtained 
by using the heat kernel method. We describe these solutions and compare them with the 
solutions of the linearized Einstein equations. We demonstrate that for D > 4 the gravi¬ 
tational potential is regular at the position of the source, and the parameter // is the scale 
of the UV regularization. We also show that for D = 4 the obtained result coincides with 
the one obtained earlier by the Fourier method and presented in the paper [7]. In section 3 
we boost the solution for a static point mass. By taking the Penrose limit with v —>■ c of 
the boosted metric, we obtain a solution of the ghost-free equations for gravitational field of 
a ’’photon” moving in D-dimensional spacetime. We show that these solutions are similar 
to the field of non-rotating gyratons [27, 28, 29]. The main difference is that a function of 
(D — 2) transverse variables, which enters the metric of the gyraton and which is a solution of 
the (D — 2) flat Laplace equation, in the ghost-free gravity becomes a solution of the (D — 2) 
operator a(— A)A. As a result, the singularity of the gyraton metric at its origin is smooth. 
In the next sections 4-5 we study the problem of the spherical collapse of null fluid in the 
ghost-free gravity. In section 4 we demonstrate how a solution for the spherical null shell 
collapse can be obtained as a superposition of a spherical distribution of gyratons, that pass 
through the fixed point and are the generators of the null shell. The method can be used in a 
spacetime with an arbitrary number of dimensions. In the present paper we illustrate it for a 
special case of the four dimensional spacetime. We use this approach to obtain a solution for 
the gravitational held of a thin null shell in the ghost-free theory. We show that a solution, 
obtained in the developed gyraton-based approach, correctly reproduces the known solution 
for the collapsing null shell in the linearized Einstein gravity. We obtain also a solution for 
the ghost-free gravity and compare it with the solution for the Einstein gravity. In particular, 
we calculate the Kretschmann curvature invariant for the solution and demonstrate that its 
singularity is smoothened. However, it remains divergent at the origin r = 0. Finally, we 
study the collapse of a spherical thick null shell and obtain its gravitational held by taking a 
specially chosen superposition of thin null shells (section 5). We demonstrate that the ghost- 
free gravitational held is everywhere regular in this model, and for small mass M the apparent 
horizon (and hence a black hole) does not form. This means that in the ghost-free gravity 
there is a mass gap for the mini black hole formation in the gravitational collapse, and, in 
this sense, its properties are somehow similar to the properties of the theory with quadratic 
in curvature corrections, discussed in [3]. Section 6 contains summary and discussions of the 
obtained results. Appendices contain details of the derivation of the thin null shell metric 
from the gyraton solutions. In this paper we use the sign convention adopted in the book 
[30]. 
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2. Non-local Newtonian gravity 

2.1 Newtonian potential in the ghost-free gravity 

In this paper we consider solutions of the linearized ghost-free gravity equations. We start 
by analyzing static solutions for a point mass in this approximation. Such a solution gives 
a modified Newtonian potential. The gravitational field of a point source in four dimensions 
was obtained earlier in [7, 15]. The authors of these papers used the Fourier method for this 
purpose. We demonstrate that this result can be obtained much easier by using the method 
of the heat kernels. This method, practically without any changes, allows one to solve a 
similar problem in the higher dimensional case. We consider a flat spacetime and denote by 
D the number of its dimensions. It is often convenient to use a number n, connected with D 
as follows D = n + 3. The static metric in the weak-field approximation in the standard D 
dimensional gravity can be written in the form (see e.g. [29]) 

ds 2 = —(1 + 2<p)dt 2 + ^1 — — (p'j dl 2 , d( 2 = 5ijdx l dxi = — ( 2 - 1 ) 

where the Newtonian potential <p satisfies the equation 

A p = 4vr Gp , G = . (2.2) 

The operator A is a standard (D — l)-dimensional Laplace operator, which in the flat (Carte¬ 
sian) coordinates is of the form 

A (...) = **(...)«■ (2.3) 

In the four-dimensional spacetime n = 1, so that G = G 4 and (2.2) takes the standard form 
of the Poisson equation for the Newtonian gravitational potential. It should be emphasized 
that there exists an ambiguity in the choice of the form of the higher dimensional coupling 
constant. We denote this constant by G D and fixed it by the requirement that the form of 
the Einstein equation is the same in any number of dimensions 

%] = 16 ^G 5 / d ° x ^ R - ( 2 - 4 ) 

We introduce another constant G, related with G D in such a way, that the form of the 

equation for the Newtonian potential p is the same for any D. For a point mass m 

p = m.5 n+2 (x ), (2.5) 


and the potential p is (see e.g. [29]) 

_ 2r(l + n/2) Gm 
^ / K n / 2 n r n 

In the four-dimensional case this expression takes the standard form p = —Gm/r. 


( 2 . 6 ) 


- 4 - 



Following [7] we write the modified ghost-free equation for the Newtonian gravitational 
potential <p, created by a point source of the mass m, in the form 

F<p = 4irGp , F = o(A) A . (2.7) 

We shall use the following form of the non-local operator o(A), proposed in [7] 

a{ A) = e- A/ ^ 2 . (2.8) 


The parameter p specifies the characteristic energy scale, where the adopted theory modifi¬ 
cation becomes important. In what follows we shall also use another parameter, s, related 
with /i as follows 


s 



(2.9) 


This is convenient, since some of the relations will contain derivatives and integrals over s, 
which do not look ’’elegant” in terms of p. On the other hand, p has more direct ’’physical 
meaning” as the mass (energy) of the UV cut-off of the modified gravity. After required 
calculations are performed one can simply substitute the parameter s in terms of p by using 
the relation (2.9). 


2.2 Heat kernel approach 

Let us discuss now equation (2.7) and its solution. We shall use the following standard 
notions. For any operator B we denote by B(x,x') its value in the x-representation 

B(x,x > ) = (x\B\x') . (2-10) 

In these notations h-function is 

6(x, x') = (x\I\x') , (2-11) 

where I is a unit operator. Using these notations one can write the equation (2.7) as follows 

J dx' F(x, x')ip(x') = AirGp(x). (2-12) 

We denote by A an inverse to F operator 

A = F- 1 = ^e sA = e sA ^~. (2.13) 

A A v ’ 

We denote 

L) = -l, K(s) = e sA , (2.14) 

so that relation (2.13) takes the form 

A = -DK(s ) = —K(s)D . 
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In the x- representation D(x, x') is nothing but a usual Green function of the Laplace operator 
and it obeys the equation 

A D = -I. (2.16) 

The other operator, K(s), is also well known. In x-representation it is just a standard heat 
kernel K n (x,x'\s) of the Laplace operator. It has the following form 

K ” (l ' I,|s)S p^75 e " X( "' l/(4S) ' <2 ' 17) 

This formula requires some explanations. We remind that we denoted by D = n + 3 the total 
number of spacetime dimensions. Since the gravitational field of a static source is static, the 
gravitational potential depends only on D — 1 = n + 2 spatial coordinates. Correspondingly, 
the heat kernel (2.17) is a function of 2 spatial points, x and x’ and A(x, x') is the square of 
the spatial distance between these points 

A(x, x') = \x — x'\ 2 . (2-18) 


The degree 1 + n/2 of expression in the denominator reflects the fact that we are working in 
a space with D — 1 = n + 2 dimensions. 

It is easy to check that the operator A can be written in terms of K (s) as follows 

/ OO 

dsK n (s). (2.19) 


Relations (2.17) and (2.19) allow one to obtain a solution of the equation (2.7)-(2.8). For a 
point mass m at point x' = 0 one has 


p{x) = 4:irmA n (x, 0) 


—4irGm / ds K n (x, 0|s) 


( 2 . 20 ) 


Before presenting explicit form of the solutions for the Newtonian potential p in the ghost- 
free gravity, let us make the following remark. The equations (2.7)—(2.8) can be identically 
rewritten in the form 

Ap = AitGp , p = K n (s)p. (2-21) 

In other words, the Newtonian potential (p in the ghost-free gravity is a solution of the stan¬ 
dard Poisson equation for the source p, which is obtained by smearing the original distribution 
p(x). For the point mass (2.5) the smeared source is 


p(x) = mK n (x,0\s). 


( 2 . 22 ) 


In this approach one may say that the ghost-free gravity regularizes the potential p by 
smearing a source, which generates it. Of course, the results obtained by the source smearing 
and by applying A operator to ”un-smeared” source p are the same. Let us present the 
results. 
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In flat three-dimensional space (D = 4, n = 1) in Cartesian coordinates 

K\(x,x'\s) = 


( 47 r.s ) 3 / 2 

4 s 

5 

1 rrf 1 

'V\ 

4ttv/A \ 

2 yf* 


Ai(x,x') = - 

v = -— erf ( T ) • 

Here erf(x) is the Gauss error function, which is defined as follows 


erf(x) = —= 
n r . 


" dt. 


(2.23) 

(2.24) 

(2.25) 

(2.26) 


Formula (2.25) correctly reproduces the expression for the gravitational potential in the lin¬ 
earized ghost-free gravity, obtained earlier [7, 15]. 

In flat four-dimensional space (D = 5, n = 2) in Cartesian coordinates 

1 


K 2 (x,x'\s) = 


A 2 (x,x') = - 


( 47 Ts ) 2 

1 

47t 2 A 


_ 

e s 


1 — e 


—V (4s) 


<P = 


Gm. 


irr- 


(1 — exp(—/r 2 r 2 /4)) . 


(2.27) 

(2.28) 
(2.29) 


There exist simple relations between K n (x, a/|s) and A n (x, x') in spacetimes with different 
number of dimensions 

\ 0 l 0 

K n + 2 (x,x'\s) = K n (x,x'\s ), A n+2 = A n (x,x '). (2.30) 

7r oa 7r a A 

Using these relations one can obtain a general expression for A n (x,x') 

n> ~ 2 - <2 ' 31) 

This relation contains a so called lower incomplete gamma function j(a, x), which is is defined 
as 


7 (a,x)= f t a 1 e t dt = x a T(a)e x ^ 
•'° k =o 


x 


T(a + k + 1 ) 


Thus in the ghost-free gravity the Newtonian field of a point mass m is 

7 (n/2,r 2 /(4s)) 


(2.32) 


(2.33) 


At small x one has 7 (a,x) ~ x a /a. Hence 

2 Gm _ „ 

y ( r = 0) — ~ n7r (n_ 1)/2 ( 4s )n /2 ' ( ' ^ 

This relation implies that the Newtonian potential of a point mass in the ghost-free gravity 
in any number of dimensions D > 4 is finite at the origin. In other words, this potential is 
properly regularized. 
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3. Non-spinning gyratons in the ghost-free gravity 

We obtain now the gravitational field of an ultra-relativistic particle in the framework of the 
ghost-free gravity. Instead of solving the modified gravitational equations for a source moving 
with the speed of light we use the following procedure which is well known in the standard 
General Relativity. We first make the Lorentz transformation of the static solution for a point 
mass m and obtain the metric for the object moving with velocity (5. After this we take a so 
called Penrose limit of the metric of the moving body. Namely, we take the limit fd —» 1, while 
keeping the energy of the object 7 m fixed. As a result one obtains an object in D-dimensional 
spacetime which was called a gyraton [27, 28, 29]. In a general case, when a static source has 
an angular momentum, the corresponding gyraton is spinning. We restrict ourselves by the 
case of non-spinning gyratons. A new element in our derivation is performing the described 
procedure not within the General Relativity, but in the ghost-free theory of gravity. 

To perform the calculations it is convenient to use the following notations for the standard 
Cartesian coordinates x^ = (t, y, £ ± ). Here t is just the time in the frame, where the source is 
at rest, y is the coordinate in the direction of motion of the source, and = (C 2 , • • •, C D 2 ) 
are the coordinates in (D — 2) dimensional plane orthogonal to the direction of the motion. 
We shall call the latter the transverse coordinates. The Newtonian potential for a point source 
of mass m in the ghost-free gravity, obtained in the previous section, takes the following form 
in these coordinates 


ds 2 

ds 2 0 

dh 2 


ds$ + dh 2 , 


n+1 


-dt 2 + dy 2 + d(l , d(l = ^(dC;) 2 , 


—2 (p 


dt 2 + - {dy 2 + d(l) 
n 


(3.1) 

(3.2) 

(3.3) 


where (p is defined by (2.33). 

To obtain a metric for a source moving with the velocity fd (along the y-axis in the positive 
direction) we make the following Lorentz transformation 


y = t(£ - fit) = | K 1 - P)v - (1 + P)u ], 

t = i(t-PQ = |[(i - P)v + (1 + P)u\. 


(3.4) 


Here (t. £, £ ± ) are Cartesian coordinates in the new inertial frame, where the source is moving. 
We also denoted 

7 = (1 - P 2 y 1/2 , u = t-£, V = t, + ^. (3.5) 

The flat metric ds^ in the new coordinates aT = {u,v,C±) is 


dsQ = —du dv + dC± • 


(3.6) 



The form of this metric remains the same in the limit /3 —>• 1, while in this limit 


y^-'yu, t~ju, X(x, 0) ~ 7 2 it 2 + Cl, (3.7) 

dP + — {dy 2 + d(l) ~ n ^ 7 2 du 2 + ... . (3.8) 

n n 

Here ... denote sub-leading in 7 terms. As a result, the metric perturbation in this limit 
takes the form 

dh 2 = Qdu 2 , <f> = ij m ( ry 2 (p ) _ (3.9) 

n 7 —►oo 

When taking this limit, we assume, as usual, that the energy of the object, 7 m is fixed (the 
Penrose limit), and we denote 

M = G'jm . (3.10) 

We also use the following relation 


lim 7 e 1 ^ = p4irsd(u). 

7—>-00 


(3.11) 


Combining these results, we finally obtain the following expression for the function 

TL 1 

$ = -87T- M5(u)A n _i (( ± ), (3.12) 

n 

where 

A.-.KJ = - [°° , W e- g/|4,) . (3.13) 

Js (47TS) 2 

Let us notice, that i(C±) is nothing, but a solution of the ghost-free gravity equations 
for a point source, reduced to D — 2 dimensional transverse plane with coordinates £ ± . For 
large /1 it reduces to a solution of the D — 2 dimensional Poisson equations. This means, that 
in this limit the obtained solution of the ghost-free gravity equations reduces to the standard 
metric of non-rotating gyratons [27, 28, 29]. 

In what follows we restrict ourselves by considering a four dimensional spacetime. For 
this case the integral, which enters the definition of has an infrared divergence at large 
s. Let us discuss this case in more detail. Let us notice that 

2 

= ye-<i /(4 ' ) = L[ Ei (l,C!/(4») + Ei(l,C!/(V)] . (3.14) 

The parameter 77 , which has the dimensionality of the length, is an infra-red cut-off parameter. 
Here Ei(o, z) is the exponential integral defined as 


Ei(a, z ) 



dxx~ a e~ xz 


The function Ei(l, z) for small z has the following expansion 


(3.15) 


Ei(l, z) = - ln(z) - 7 + ... . 


(3.16) 
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Here 7 = 0.5772156649 is the Euler constant, and ... denote the terms vanishing in the z = 0 
limit. Assuming that i] 2 is large and using (3.16) one can write 


A 0 = ^ [ln(Cl/r? 2 ) + 7 + Ei(l, C!/(4s)] . (3.17) 

The corresponding expression for the metric is 

ds 2 = —dudv + d£± + dh 2 , dh 2 = &du 2 , (3.18) 

$ = -4 M5(u)F((l) , (3.19) 

HCl) = HCI/V 2 ) + 7 + Ei(l, C?/(4s)) • (3.20) 

In the limit fx —> 00 , when the ghost-free gravity reduces to the Einstein theory, 

F(C?) = ln(Ci/r? 2 )+7- (3-21) 


This result reproduces the well known Aichelburg-Sexl solution [31] (see also [32]). 

Let us notice that the solution (3.20) contains an arbitrary parameter, the infrared cut¬ 
off parameter ij. However, the change of this parameter can be easily absorbed into the 
redefinition of the advanced time v. Hence, this ambiguity just reflects freedom in the gauge 
choice. In order to demonstrate this, let us consider a metric of the form 


ds 2 = —dudv + d £ 2 + f(u)du 2 . 

(3.22) 

Let us define a new coordinate 


v = v — / duf{u). 

(3.23) 

Then one gets 


ds 2 = —dudv + . 

(3.24) 


This confirms our above conclusion, that the ambiguity in the choice of the cut-off parameter 
r] can always been absorbed in the change of the coordinates. In what follows we use this 
option and simply put i] 2 = 4s. For this choice 


F(*) = ln(z) + 7 + Ei(l,z), z = C 2 J(As). (3.25) 

For this choice the expansion of the function F(z) for small z takes a very simple form 

F(z) = z-^z 2 + 0(z 3 ). (3.26) 


4. Null shell collapse 

4.1 Null shell as a superposition of null gyratons 

We use now the above described gyraton metric in order to study the gravitational collapse 
in the ghost-free gravity. Namely, we consider a collapse of a spherical null shell. As earlier 
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we use linearized gravitational equations. For simplicity, we restrict ourselves by 4D case. 
Instead of solving the corresponding equations we shall use the following trick. Let us notice 
that a sum of solutions of the linearized theory is again a solution. In the linear approximation 
the solution has the form 

ds 2 = dsQ + dh 2 . (4.1) 

Here ds^ is the flat metric and dh 2 is a perturbation. 

Let us consider a set of gyratons passing through a chosen point O of the spacetime. 
We use the Cartesian coordinates (t, X,Y, Z) in the the Minkowski spacetime, and identify a 
point O with the origin of the coordinate system O = (0,0, 0,0). We denote by (e x , e Yl e z ) 
unit vectors in the directions of the axes X, Y, and Z, respectively, and by n a unit vector 
in 3D space in the direction of the motion of a fixed gyraton. One has 

n = sin a (cos j3 e x + sin f3 e Y ) + cos a e z . (4-2) 

Here (a, f3) are standard coordinates on a unit 2D sphere. Let a be an index enumerating the 
gyratons and the metric perturbation created by such a gyraton is dh 2 . Then the perturbation 
(dh 2 ) created by a set of the gyratons has the form 

(dh 2 ) = Y J dh 2 a . (4.3) 

a 

We assume that all gyratons have the same energy. We also take a continuous limit of 
the discrete distribution of the gyratons and assume that such a distribution is spherically 
symmetric. Thus, we write 

^ /*7T /*27T 

(dh 2 ) = — j da sin a J d/d dh 2 a ^ . (4.4) 

An extra factor (4vr ) _1 reflects that we use averaging over a unit sphere. 

As a result of the averaging, the source of the metric dh 2 is a thin null shell located at 
the null cones r± 

t = ±^X 2 + Y 2 + Z 2 . (4.5) 

For the sign minus, is a null cone with apex at O, which describes a collapsing spherical 
null shell. For the sign plus, the null cone describes an expanding null shell. Our starting 
point is gyraton metric (3.18), which we rewrite in the form 

dsl = -dt 2 + df + dC 2 ± , dh 2 = <f>{dt - d £) 2 , $ = -4 MF((l)5(t - £). (4.6) 


Here 


F = 


In (Cl/ 7 ? 2 ) j for the Einstein theory; 

Ei(l, Cl/4s) + 7 + ln(£^/4s), for the ghost-free gravity. 


(4.7) 


Let us notice that this is a metric of a gyraton without spin. Namely these gyratons will be 
considered in this paper. We also found convenient to use the following terminology. We call 
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a worldline of a gyraton a null string. Its projection to the 3D space, the gyraton trajectory, 
is an oriented straight line. A unit orientation vector n along this line determines a direction 
of the motion. Quantity £ is the coordinate along the trajectory in the direction of motion 
of the ’’photon”, and £x = (C 11 C 2 ) are Cartesian coordinates in the 2D plane orthogonal to 
this direction. We study a spherically symmetric distribution of the gyratons, which has the 
property that the null strings, representing them, intersect at a single spacetime point O. We 
also choose the parameter £ along each of the strings to vanish at this point O. Consider a 
point P = ( X , Y, Z) of the 3D space. There exist exactly two gyraton trajectories, passing 
through this point. Let (a + ,/3 + ) be angles of a unit vector n + (see (4.2)) along the line, 
connecting the origin with the point P. The parameter at P for such a trajectory is 
positive. The direction vector for the second trajectory is n_ = — n + and its angles are 
(a- = ir — a + , (3- = 7T + /3 + ), while the corresponding coordinate £_ = — £+ is negative. 

It is convenient to perform the calculations of ( dh 2 ) in two steps. First we introduce the 
following objects 

T^\y ± ) = 5 2 (y ± - { ± )5(t - $(dt - dtf , (4.8) 

1 /*7T 

{T(y±)) = — j da sin aJ d/3 T (a>/3) (y x ). (4.9) 

For y± = 0 

r (a,/3) = T (*,P) dx » dx v (4.10) 

is the stress-energy tensor of a gyraton moving in the direction n (see (4.2)). Similarly, for 

y± = 0 

M(T) =T„ v dx lx dx v , (4.11) 

where T /IU is the stress-energy tensor of a spherical null shell, constructed from gyraton null 
strings. Next, we use (T(y ± )) to find the metric perturbation (dh 2 ) for the thin null shell of 
mass M 

(dh 2 ) = M J dy ± F(y 2 ± )(T(y ± )) . (4.12) 

It should be emphasized that the quantities dh 2 Q m and T^ a, d)(y ± ), which enter (4.4) and 
(4.9), must be first written in the coordinate system, that does not depend on the particular 
value of the parameters (a,/3). We use the Cartesian coordinates (X. Y. Z) for this purpose. 
Thus we need first to establish relations between the gyraton associated coordinates (£,Cx) 
and the Cartesian coordinates (A, Y, Z ). This problem is solved in the appendix A. After this 
we need to calculate the integral over the sphere, which enters relation (4.9) for the average 
value (T(y ± )). The details of these calculations are collected in appendix B. 

4.2 Stress-energy tensor 

The relation (B.18) can be used to find the stress-energy tensor of the null shell constructed 
from null strings representing a set of gyratons. Using (4.9) and (B.18) one gets 

(T(y±)) = 4 1 * [sma + 5(t-t + )(dt-d{, + ) 2 + sma-8(t-t-)(dt-d£-) 2 ] y± . (4.13) 
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Taking the limit y ± —> 0 in this relation and using relations (4.11) and (B.14) one obtains 

M(T) = T^dx^dx" = ~^[5(u)du 2 + 5(v)dv 2 ), (4.14) 

where u = t — r, v = t + r. This relation correctly reproduces the expected expression for 
the stress-energy tensor of the null shell. It is a superposition of the stress-energy tensors of 
contracting and expanding spherical null shells of mass M. One can easily solve the linearized 
Einstein equations for such a null-shell problem. The solution is well known and simple. Inside 
both the collapsing and expanding shells the spacetime is flat, while outside them the metric 
is a linearized version of the Schwarzschild metric 

(dh 2 ) = ^-(dt 2 + dr 2 ). (4-15) 

r 

Our next goal is to reproduce this result by using the representation (4.12). This will provide 
us with a useful test of the validity of our approach. 

4.3 Averaged metric 

We have all the required expressions to perform the calculations. However, one can greatly 
simplify the problem using the following observation. Since the distribution of the null strings 
representing gyratons is spherically symmetric, the corresponding averaged metric (dh 2 ) must 
also have this property. Hence, it can be written in the form 

(dh 2 ) = httdt 2 + 2 htrdtdr + h rr dr 2 + Hdu 2 . (4-16) 

Here the metric coefficients htt , ht r , h rr and H are functions of t and r, and du 2 is the metric 
on a unit sphere. The spherical coordinates (r, 9, </>) are related with the Cartesian coordinates 
(A, Y, Z) as follows 

X = r sin 9 cos (/>, Y = r sin 9 sin <f >, Z = rcos9. (4-17) 

In order to find the metric (4.16) it is sufficient to obtain its value near a single spatial point. 
We choose such a point as follows: P x = (X = r, 0,0), where X > 0. For this choice the 
metric (4.16) reduces to 

(dh 2 ) = h tt dt 2 + 2 h tr dtdX + h rr dX 2 + — 2 ( dY 2 + dZ 2 ) . (4.18) 

We perform now calculations of (dh 2 ) at Px and by comparing it with (4.18) we find the 
expression for the metric perturbation (4.16). One has near P x 


Q± = ±^x 2 yj , e± = ±V^ 2 - vi, 

(4.19) 

■ Vk 

sin a± = —— , cos a± = — —— , 

Q± ' Q± 

(4.20) 

sin/3± = ^, cos/3± = ^, 

(4.21) 

d£± = JtdX + -±-(t±?£dY-y k dz) . 

(4.22) 
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In order to obtain the metric perturbation ( dh•?) one must to calculate the integral over 
(y p ,y k ) in (4.12). We introduce polar coordinates in the (y p ,yk)~ plane: 


We also write 


y p = p cos ip , Vk = P sin ip 


(T(y ± )) = T + + T~ , 


T± = 


AttQ- 


- $(t ~ €±) (dt - d£±)‘ 


(4.23) 

(4.24) 

(4.25) 


Because of the presence of 4-function the term T + is non-zero only for t > 0, while the other 
term T~ is non-zero for t < 0. We write expression (4.12) in the form 


One has 


(dh 2 ) = -AM / dppF(p 2 )(T + + T ~), 

Jo 

T ± = dip . 

Jo 


(dt - dtp±) 2 = A± + ^ + ... , Q 2 = X 2 - p 2 cos 2 ip . 


Q 

A± = (dt — ^dX) 2 , B± = (a± cos ip + bsinip) 2 , 
si. 


a± = 


t ±pdY 
X 


, b = —pdZ , 


(4.26) 

(4.27) 

(4.28) 

(4.29) 

(4.30) 


and dots indicate terms linear in y p and ijk, which do not contribute to the integral over ip 
(4.27). Thus one has 


t* = S -^{a±k 1 + k 2 ), Jf, = ["%, a - 2 = 


/»27T 


dip B± 

~Q r ~ 


Integrals K\ and K -2 can be easily taken with the following result 

2t r 

xjW^J 2 ' 

7r[(a 2 ± + b 2 )X 2 - b 2 p 2 } 


K i = 
K-2 = 


irp 


X 3 (X 2 - p 2 ) 3 / 2 X 3 y/X 2 - p 2 
4-function in (4.31) can be written as 


(dY 2 + dZ 2 ). 


(4.31) 

(4.32) 

(4.33) 


S(t v 7 X 2 - P 2 ) = —S(p - Vx 2 - t 2 ). 

p 


This relation also implies that 


t = ±y/X 2 - p 2 


(4.34) 

(4.35) 
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After integration in (4.26) the terms T± give similar contributions, so that finally one obtains 
the following result 


{dh2) = jMgMj 

X 


dt-^dX) + 


X 2 - t 2 
2X 2 


( dY 2 + dZ 2 ) 


(4.36) 


Comparing (4.36) with (4.18) we finally get 




t \ r 2 — t 2 2 

dt - dr H- du 2 

r 2 


(4.37) 


This expression is valid everywhere in the domain r > |t| for both positive and negative time 

t. 

Using GRTensor program one can check that the Ricci tensor for this perturbation of the 
flat metric in the linear in M approximation vanishes everywhere outside the null shell. This 
provides one with a good test of the correctness of the performed calculations. 


4.4 A case of linearized Einstein equations 

For the linearized Einstein equations 

F(r 2 - t 2 ) = In ^ . (4.38) 

The corresponding expression (4.37) looks quite different from the expected answer (4.15). 
However one can use the gauge freedom 


to transform (4.37) into the expected form. It is sufficient to choose 


(4.39) 


^ = (^,^, 0 , 0 ), 

2Mt. 

Vt = -(1 + t/r ) In 


r 2 — t 2 


7j- 

M , o o. , ( r 2 — t 2 


r — t 2 Mt 

+ 2 M In-1- 

r +1 r 


V r = -72 ( r ) ln 


7J~ 


After this gauge transformation one gets 


(dh 2 ) = ^^-(dt 2 + dr 2 ). 
r 


(4.40) 

(4.41) 

(4.42) 

(4.43) 


Let us notice that the infrared cut-off r/ 2 in (4.38) does not enter the final result and the 
change of this parameter is simply absorbed into the redefinition of the gauge held V fl . Let 
us mention also that the Kretschmann invariant 


n 2 = 


(4.44) 



in its lowest non-vanishing order is 

, 48 M 2 

77 2 = —— 

as it must be for the Schwarzschild metric. Let us emphasize that the same result can be 
obtained directly from the perturbation of the metric in the form (4.37). 



4.5 Ghost-free case 


One can easily repeat similar calculations for the an arbitrary function F(z), where where z = 
r 2 —t 2 . We used the GRTensor program for this purpose. The calculations are straightforward. 
However, the intermediate formulas are quite long. That is why we do not reproduce them 
here. Let us only present the expression for the Kretschmann invariant in its lowest order for 
such the metric in the ghost-free gravity 

77 2 = - Q{z) , Q(z) = 2F' 2 z 2 + 2FF'z + T 2 , 7 = F’ . (4.46) 

For the ghost-free theory with a(D) = exp(— \3/fi 2 ) one has 


F(z) = Ei( 1, z/(4s)) + 7 + In (z/ (4s)), 

F{z) = \ (l - e"^ 4 *) , 

Q(z) = ^7 (W + 4sz + z 2 )e~ z/2s - 4s(4s + z)e~ z / 4s + 8 s 2 ) . 


Using the expansion (3.26) of F for small z 


z z 

F ^ = 4 ^ “ 64 ? 


+ • 


one finds 


2 


77 2 = 


3 M 2 z 
s 2 r 6 


2«2 


+ . . . = 


3 M 2 P 

S 2 ?’ 4 


+ 


(4.47) 

(4.48) 

(4.49) 

(4.50) 

(4.51) 
< r, is 


Here f3 2 = 1 — t 2 /r 2 is a dimensionless parameter which outside the shell, where 
less or equal to 1. This means that the curvature for the linearized ghost-free theory in a 
case of the collapse of the spherical null shell is weaker than the singularity for the linearized 
Einstein equations. However, the ghost-free gravity solution still remains singular at least for 
the chosen scheme of the regularization. 


5. Spherical ’’thick” null shell collapse: Results 
5.1 ’’Thick” null shell model 

In the previous section we discussed the gravitational field of a spherical null thin shell. This 
shell represents a spherical 5-type distribution of the energy. The mass of the shell is M. It 
collapses with the speed of light and shrinks to zero radius at the moment of time 7 = 0. In the 
linearized theory the corresponding perturbation of the background flat metric is ( dh 2 )(t,r ). 
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Figure 1: Thick null shell collapse 

Certainly, such a model is an idealization. To study more realistic model of the collapse we 
assume that the spherical collapsing null fluid is represented by a pulse, which is not infinitely 
sharp in time, but has final time duration. We characterize its profile by a function q{t). The 
meaning of this function is the mass density per a unit time, dM/dt, arriving to the center 
r = 0 at time t. The total mass of such ’’thick” shell is 

M = j dt q{t ). (5-1) 

Using the linearity of the equations one can write the corresponding solution (( dh 2 )){t,r ) as 
the superposition of ( dh 2 ){t,r ) perturbations as follows 

((dh 2 )){t,r) = J dt'q{t')(dh 2 ){t — t',r). (5-2) 

To simplify the calculations we further specify the model. Namely, we choose q{t) to 
be a step function, which has a constant value M/b during the interval t E (—6/2, 6/2), and 
which vanishes outside this interval. For such a step function a solution has different form 
in the different spacetime domains (see Figure 1). Let us describe first these domains. It is 
convenient to use the advanced time, v = t + r, and the retarded time, u = t — r, coordinates. 
In the domains T_, where v < —b/2, and T + , where u > 6/2, the spacetime is empty, and 
the metric is flat there. In the domain iV_, where v E {—b/2, 6/2) and u < —b/2, one has 
only in-falling null fluid flux, while in N + , where u E {—b/2, b/2) and v > b/2 one has only 
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out-going null fluid flux. In the domain, /, where v G (—6/2, 6/2) and u G (—6/2, 6/2) one has 
a superposition of the in-coming and out-going null fluid fluxes. And finally, in the domain 
/?, where v > 6/2 and u < — 6 / 2 , the spacetime is empty. 


5.2 Gravitational field in the /-domain 


Let us consider first /-domain. Denote by (t,r) the coordinates of the point of ” observation” 
in this domain. A thin null shell, crossing the center r = 0 at time t', contributes to the 
integral (5.2) only if i! G (t — r, t + r). Let us denote x = t' — t, then the formula (5.2) takes 
the form 


9 M r 

(( dh 2 )) (t, r) = —— J dxH (r, x ), 

H(r,x) = F(r 2 — x 2 ) f/f 2 — 2-drdt+ ^ dr 2 + -(?’ 2 

\ r r z 2 

Thus one has 

J 0 {dt 2 + ^r 2 /o; 2 ) + ./ 2 ( rfr 2 - ^w 2 ) 
2 r z 2 


(( dh 2 ))(t,r ) = - 


2M 

hr 



(5.3) 

(5.4) 


(5.5) 


Here 


/ r 

dx x n F{r 2 — x 2 ). 

-r 

The linear in x term in (5.37) disappears as a result of the integration. 


(5.6) 


5.2.1 Linearized Einstein theory 

For this case 


Calculations give 


F(z) = In z + C , C = — In f ? 2 , 

Jo = 4r(ln2 + lnr — 1 ) + 2Cr , 

4 2 

J 2 = -r 3 (31n2 + 3 ln r — 4) — Cr 3 . 
9 3 


252A/ 2 
77 2 =-, 


M = M/b. 


5.2.2 Ghost-free gravity case 

One has 


F(z) = Ei(l, z/4s) + 7 + ln(z/4s). 
The Tailor expansion of this function for small 2 is 


z z 

F ( : ) = 4 ^ “ 64 ? 


+ 


1152s 3 


+ ... • 


(5.7) 

(5.8) 

(5.9) 

(5.10) 


(5.11) 


(5.12) 
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Using this expansion one gets for small r in the /-domain 

T r 3 (420s 2 — 21sr 2 + r 4 ) 

0 “ 1260s 3 + “' ’ 

_ r 5 (756s 2 - 27sr 2 + r 4 ) 

2 “ 11340s 3 + 

Calculations of the Kretschmann invariant 1Z 2 for small r in the /-domain give 


77 2 = 


32 M 2 32 M 2 r 2 


3s 2 


9s 3 


+ ... . 


(5.13) 

(5.14) 

(5.15) 


Here, as earlier, M = M/b. This relation shows that for the thick shell model the curvature 
remains finite at r = 0. 

5.2.3 (Vr) 2 invariant 

There exists another useful invariant for the spherically symmetric geometry. This invariant 


is 


1 


(V r ) = , f = gee- 


(5.16) 


A line in the (t. r ) plane, where this invariant vanishes, is an apparent horizon. Using GRTen- 
sor one find that in the leading order in M this invariant in the /-domain is 


(Vr) 2 = 1 - 


2 Mr 2 
sb 


+ 


Since in this domain r < b, one has 


(Vr) 2 > 1 - 


2 Mb 


(5.17) 


(5.18) 


This relation means that for given s, which is the square of the UV cut-off parameter ji of 
the ghost-free theory, s = g~ 2 , and fixed duration b of the pulse the mini-black hole is not 
formed if the mass M is small enough. 

5.3 Gravitational field in the R domain 

The perturbation of the gravitational field in the R domain is described by the following 
expression 


« dh 2 ))(t,r ) = - 2 ^ r 


1 


2 dt dr 


Jo(dt + -r 2 du ) - Ji 

2 r 


,dr 


1 


+ ^ 2 (T2- 


where now 


Here 


f‘t-\-b /2 

J n = dx x n F(r 2 — x 2 ) 

Jt—b/2 


F(z ) = F 0 (z) + A F(z ), F 0 (z) = 7 + ln(z/4s ), 
A F(z) = Ei(l, z/4s). 


(5.19) 

(5.20) 


(5.21) 

(5.22) 
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We denote by <7° and A J n the contribution to J n of the terms Fq and A F, respectively, so 
that one has 

Jn = Jn + A J n . (5.23) 

We consider the perturbed metric at large r in the R domain. We also assume that 
r 3> \t ± b/2\ and use the following asymptotic form of the function A F(z) for large z 

A F(z) = ^V z / 4s (l + 0{l/z)) . (5.24) 

Thus one can use the following approximation 

A F(r 2 — x 2 ) « f(r)k(x) , f(r) = -^e~ r2 ^ 4s ' > , k(x) = 4 se x2A4A . (5.25) 

Using this approximation and (5.20) one obtains 

A J n = f(r)A n (t), A n (t) = (4s)^ f dyy n e y2 , C± = t ± b / 2 . (5.26) 

Jc- Z V- S 


Calculation of is straightforward and <7° can be easily obtained by using the Maple pro¬ 
gram. One also obtains the following expressions for A n (t ) 


Ao{t.) 

A i(t) 

A 3 (t) 


= — (erf(ir + ) + erf(ir_)) , 

= 8s 2 (e T + - e T ~ 

= 4 s 2 \2i\fns (erf(u+) + erf(ir_)) + (2 1 + 6)e T + 


(2 1 - b)e T - 


(5.27) 

(5.28) 

(5.29) 


Here 


±2 t + b 


^ 4 yfg ’ 

and erf is the Error function, which is defined as 

2 f x 9 

erf(x) = —= / dtexp (—t ). 
V 71 Jo 


(5.30) 


(5.31) 


Let us emphasize that in spite of the presence of the imaginary unit i in the above formulas, 
expressions for A n (t) are real, as they should be. After straightforward calculations by using 
the GRTensor program, we obtain the following expression for the Kretschmann tensor in the 
leading M 2 order 


o 48 A7 2 A o 
77 2 = —g-h A77 2 


2 2A7 2 

A7Z ~ ~ s^Io e 

OA A ^ I A \~,4 


/( 4s ) w , 

W = — 3Ao?’ 6 + (4Aos 2 + 8Ais — 34Aos + A 2 )r 4 
+(80Ais 2 - 12A 2 s 2 - 160A 0 s 2 + 14A 2 s)r 2 + 56A 2 s 2 . 


(5.32) 

(5.33) 


Keeping the first term in W, which contains the highest power of r, we obtain 

A7Z 2 ~ e _r 2 /(4s) ( er f(j r ) + erf(*r_)) . 

V-sor 4 


(5.34) 
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Once again, in spite of the presence of i the answer for A1Z 2 is real. 

Expressions (5.32) and (5.34) imply that in the R domain for a fixed value of time t and 
r —>• oo the curvature exponentially fast reaches its asymptotic value, which is equal to the 
Schwarzschild curvature of a spherical object of mass M. 


5.4 Gravitational field in the A/t-domains 

Let us briefly discuss now the gravitational held in the N± domains. We focus on the incoming 
1V_ domains. The case of the N + is similar. Let (t, r) be coordinates of the ” observation” 
point. We denote v = t + r and u = t — r. The condition that the ’’observation” point is in 
the 1V_ domain reads 

v G (—6/2,6/2), u G (—oo, —6/2). (5.35) 

The position t' of the apex of the null-shells, that give non-zero contribution at the point 
of the ’’observation”, obeys the conditions t' G (—6/2, t + r). As earlier, we denote x = t — t!. 
Then 


„ 2 M [ t+b/2 

{(dh 2 }) (t, r) = —— J dxH(r, x ), 

H(r, x) = F(r 2 — x 2 ) f (dt — -dr) 2 + -( r 2 — x 2 )duj “ 
V r 2 


Thus one has 


(( dh 2 ))(t,r ) = - 2 ^ r 


Here 


Jo (dr + - r z dcu z ) - Ji 

z r 


J n = dx x n F(r 2 — x 2 ). 

J—r 

In the linearized Einstein gravity one uses the following expression for F 

F(r 2 — t 2 ) = In 


2 dt dr .dr 2 1 2 

+ J 2(W - n duJ 


r 2 2 


r 2 — t 2 


rj* 


(5.36) 

(5.37) 

(5.38) 

(5.39) 

(5.40) 


The required integrals (5.20) can be easily calculated. We do not reproduce here these results 
since the obtained expressions are quite long. We present here only final expression for the 
Kretschmann tensor, which we obtained by using th GRTensor program 

77 2 = l,V ^" >(r) , M(v) = M(v + 6/2). (5.41) 

This is an expected result. Indeed, before the collapsing null fluid meets the outgoing flux, 
that is in the 1V_-domain, one has simply the Vaidya solution with mass M(y) = ^-(v + b/2), 
where v G (—6/2, 6/2). For a chosen model M{y ) linearly grows from zero (at v = —b/2) till 
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its maximal value M. Calculations for the metric perturbation in the N + are similar and give 
the following result 

n 2 = lN '^ (A) , M(u) = M(u + 6/2). (5.42) 

In the case of the ghost-free gravity the curvature is modified in the narrow ’’skin” 
domains close to v = ±6/2 (for 1V_ domain) and close to u = ±6/2 (for N + domain), while 
outside of them and at large r the ghost-free theory corrections only slightly modify the above 
described Vaidya metric. 

6. Summary and discussion 

Let us summarize and discuss the obtained results. We remind that we study the modification 
of the solutions of the Einstein theory in the framework of the ghost-free theory, proposed in 
[7, 8, 10, 11]. We focus mainly on a special type of such a theory in which a free field operator 
□ is replaced by the non-local operator of the form exp(—D/// 2 ) □. We study solutions of 
such a theory in the linearized approximation. We first study a gravitational field of a point 
mass in this theory. Such a problem was solved in the four dimensional space time in [7, 15] 
by using the Fourier methods. We demonstrate that the propagator of the linearized ghost- 
free theory in any number of spacetime dimensions is directly related to the heat kernel in 
this space. Using the heat kernel approach we obtain solutions for the gravitational field of 
a point mass in D-dimensional spacetime. We demonstrate that these solutions are always 
regular at the position of the source. In 4D spacetime the obtained solution coincides with 
the result presented in [7, 15]. 

In the second part of the paper we study solutions of the ghost-free gravity for ultra- 
relativistic sources. We first obtain ghost-free analogues of gyraton solutions [27, 28, 29] 
in the ghost-free gravity. For this purpose we boost the obtained ghost-free solution for 
a point mass, and by taking the Penrose limit, we found a solution for a D-dimensional 
spinless “photon”. We demonstrate that the corresponding metric is similar to the metric 
of a spinless gyraton. The main difference is that a solution of (D — 2)-dimensional flat 
Laplace equation for a point charge, that enters the gyraton metric, is modified. Namely, the 
function of the transverse variables becomes a solution of (D — 2) ghost-free modification of 
the corresponding Laplace operator. Using this result we obtain spinless gyraton solutions 
of the ghost-free gravity in D-dimensional spacetime and demonstrate that their transverse 
singularities are regularized. 

Finally, we study a spherical gravitational collapse of null fluid in the framework of the 
ghost-free gravity. As earlier, we restrict ourselves by working in the linearized approxima¬ 
tion. Since a solution of the ghost-free equations for the dynamical problem seems to be 
complicated, we used the following trick based on the results, obtained earlier in this paper. 
We used a fact that in the linearized theory a superposition of any solutions is again a so¬ 
lution. To obtain the gravitational field for a collapsing spherical thin null, we ’’construct” 
it as a superposition of the gyraton metrics for a spherically symmetric distribution of the 
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gyraton sources moving along a null cone, representing the shell. Such gyratons intersect at 
a single point of the spacetime, which is an apex of the null cone. In the adopted linearized 
approximation the gyratons, forming the thin shell, cross this point simultaneously without 
interaction, so that after passing the apex point they form an expanding spherical null shell. 
The gravitational field for such null shells is obtained by averaging of the single gyraton 
metric over homogeneous spherical distribution of the gyratons. We demonstrated that by 
using this approach one correctly reproduces the stress-energy tensor for both, collapsing and 
expanding thin null shells. We checked that the obtained solution in the linearized Einstein 
gravity correctly reproduces the expected result. Namely, the gravitational field inside both, 
collapsing and expanding shells vanishes, while outside of them it is time independent. The 
latter field is nothing but a linearized version of the Schwarzschild metric with mass M. Af¬ 
ter this, we obtained a similar solution for the ghost-free gravity. To study its properties we 
calculated the Kretschmann scalar for it and demonstrated that its singularity at r = 0 is 
smoothened, but still is present. 

The model of an infinitely thin null shell is certainly an idealization. One can expect 
that in a consistent theory the non-locality cannot be a property of the gravitational field 
only, but similar non-locality must be present in the description of matter and other fields. 
For this reason a “physical” collapsing shell must have natural finite thickness, which may 
regularize the curvature singularity of the thin shell. In order to check this assumption, we 
study a case of the thick null shell collapse. In the linearized theory a corresponding solution 
is obtained by a superposition of the thin-shell solutions, that is by averaging these solutions 
with different positions of their apexes at r = 0. As a result, we construct a solution for a 
thick null shell with an arbitrary distribution of the collapsing mass M = M ( v ) at the spatial 
infinity. We considered in detail a special model, where the mass M ( v ) is a linear function of 
the advanced time v = t + r, such that M vanishes for the moments v < —6/2 and v > 6/2, 
and remains constant inside this interval. We demonstrated, that in the linearized Einstein 
gravity the solution is a superposition of two linearized Vaidya metrics, for the incoming 
and outgoing null fluid fluxes. The corresponding Kretschmann scalar is 1Z 2 = 25 ^ , where 
M = M/6. It is singular at r = 0. 

In the linearized ghost-free gravity the corresponding solution is modified. The main 
new feature is, that the metric is regular at r = 0. As a result, its Kretschmann scalar 
1Z 2 = ’/M' is finite at r = 0. We also calculated the invariant (Vr) 2 and demonstrated for 
the collapse of a thick null shell of the small mass M the apparent horizon is not formed. 
This result can be interpreted as follows: in the ghost-free gravity there exists a mass gap 
for the black hole formation in the gravitational collapse. This result is similar to the result 
obtained in [3] for the null shell collapse in the theory gravity with quadratic in the curvature 
corrections. Recently it was demonstrated that the mass gap for mini black hole formation 
is a common property not only of the ghost-free gravity, but also of a wide class of higher 
derivative theories of gravity [33]. These theories contain a mass scale parameter /i, which 
plays the role of the ultra-violet cut-off. As a result, if the mass M of a collapsing object 
obeys the relation M/j < 1, an apparent horizon is not formed. The presence of such a scale 
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parameter, differs these theories from the classical Einstein gravity, where the mass gap is 
absent and the black-holes of arbitrary small mass can be formed. A well known consequence 
of this is Choptuik [34] type universal scaling properties of near-critical solutions. 

Let us make a few general remarks concerning solutions of the ghost-free gravity equa¬ 
tions. The heat kernels, which we used to construct solutions for a static source, are well 
defined in the space with the Euclidean metric. The reason is that the Laplace operator is 
non-positive definite. However this property is not valid for the box operator. In this paper 
we used the method based on the gyraton solutions to overcome this problem for a special 
case of the null shell collapse. It would be interesting to investigate solutions of the linearized 
ghost-free gravity with the form-factor a = exp(—□ //i 2 ) for arbitrary moving sources, say, 
for the emission of the gravitational waves in such a theory, and the back-reaction of this 
radiation on the accelerated objects. Another proposed option is to modify this form-factor, 
and instead of □ to consider higher powers of this operator in the exponent. A simplest 
example is a = exp(o0 2 ). It will be interesting to study such modifications. 

In the present paper we study solutions of the linearized equations of the ghost free 
theory. Let us emphasize, that this is sufficient for demonstration of the existence of the 
mass gap for mini black hole formation in the adopted model. The reason is that if such a 
solution for small mass M is regular and the corresponding perturbation of the flat metric is 
uniformly small, then one can expect that the higher order corrections to this solution can be 
neglected. A natural question is what happens in the collapse of large mass M, that is when 
M/i > 1. It would be highly interesting to analyze solutions of the ghost-free gravity in this 
regime [15, 16]. It is natural to assume that curvature remains finite inside the black holes 
and the limiting curvature conjecture is satisfied [35, 36, 37]. These was a lot of discussions of 
non-singular models of black holes. One of the option is that the apparent horizon is closed 
[3] (see also [38, 39, 40] and references therein). Another option is a new universe formation 
inside a black hole [41, 42], This option was also widely discussed in the literature. One may 
hope that proposed ghost free modifications of the Einstein gravity, which are ultraviolet 
complete and asymptotically free, would allow one to answer intriguing questions concerting 
the structure of the black hole interior. 
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Figure 2: Relation between the Cartesian and the gyraton based coordinates 

A. Useful geometric relations 

In order to find this coordinate transformation let us introduce a couple of new useful unit 
vectors. The first one is a vector k which is orthogonal to n and is directed from the gyraton 
trajectory to the Z-axis (see Figure 2). It can be written as a linear combination 


k = k z e z + k n n . 

Using relation (n, e z ) = cos a one finds 


k = 


1 


sm a 


-(e z — cos an). 


(A.l) 


(A.2) 


(A.3) 


The second unit vector p is orthogonal to both n and k. It can be written as 

p = [n x k] = - [n x e z ] . 

sin a 

Hence this vector is orthogonal to e z and lies in the X — T-plane. It is easy to check that 

p = sin /3 e x — cos (3 e Y ■ (A.4) 

Consider a point P, which has the Cartesian coordinates (X, Y, Z), so that the vector P 
connecting the origin point O with P is 


P — X&x + I e Y + Ze z . 


(A.5) 
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The same vector can be written in the gyraton’s frame as follows (£ x = (CCp)) 


P = £n + ( k k + C P P • (A.6) 

By comparing (A.5) and (A.6) and using relations (4.2), (A.2) and (A.4) one finds 


X = cos /3 (£ sin a — Ck cos a) + Cp sin /3 , 

Y = sin f3 (£ sin a — Ck cos a ) — Cp cos /3 , (A.7) 

Z = C cos a + Ck sin a . 

We shall also use the following inverse relations 

C = sin a (cos f3X + sin f3Y) + cos a Z , (A.8) 

Ck = ~ cos a (cos /3X + sin j3Y ) + sin a Z , (A.9) 

Cp = sin/3 X — cos f3Y . (A. 10) 

Using these results it is easy to show that 

dC 2 + d(l = dX 2 + dY 2 + dZ 2 . (A.11) 


B. Evaluating integrals 


We assume that Cp takes the value y p and use relation (A. 10) to find a corresponding value 
of the angle (5. By solving this equation one gets 


sin (/?±) = 
cos (/3±) = 


y p X + YQ± 
X 2 + Y 2 ’ 
-y v Y + XQ± 
X 2 + Y 2 


(B.l) 
(B.2) 
(B.3) 

These relations show that for each value of y p there exist two different values of the angle /?, 
which we denoted by f3±. Using (B.l) it is easy to check that 


Q ± = ±JX 2 +Y 2 -y 2 


cos f3±X + sin j3±Y = Q± , 

and the relations (A.8) and (A.9) take the form 

C = sin aQ± + cos aZ , 

Ck = — cos aQ± + sin aZ . 

Let us calculate the following integral 


[• • -]y P ~ 2 tt f 0 Cp) (•■•)■ 


(B.4) 


(B.5) 

(B.6) 


(B.7) 
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Using the relation (A. 10) one gets 

d( P = \dp(p\d/3 = Q d/3 , 


Q = Q+. 


(B.8) 


Using this relation one rewrites (B.7) as an integral over C, p which can be easily taken with 
the result 


\Vv 


1 

2irQ 


(B.9) 


■)y p + (• • Oj/pJ 

The expression in the right-hand side is a sum of the contributions for two values of /3± 
corresponding to the same y p . The quantity Cp hr each of these terms should be substituted 
by y P - We use a similar trick as earlier to calculate an integral over the angle variable a. 
First, we use the relation (A.9) with C,k = Uk to find the angle a. One gets 

VkZ + Q+\/r 2 -yi 


sina-t = 


r 2 - y 2 


cos ay- = 


~VkQ± ± Z^Jr 2 - y\ 


r 2 - y p 
r 2 = X 2 + Y 2 + Z 2 , 


2 2,2 
y± — y P + 2 /?; • 


One also has 


£± = ±y/r 2 - 


yl 


(B.10) 

(B.ll) 

(B.12) 

(B.13) 


In order to obtain this results one needs to solve quadratic equations, which have two 
solutions. We single out a solution as follows. A line yk = y p = 0 coincides with the gyraton’s 
trajectory and a and /3 are spherical angles of its direction. For y^ = y p = 0 relations (B.10), 
(B.ll) and (B.13) give 

yj X 2 + Y 2 Z 

sina-t = -, cos ay- = ± — , = ±r . (B.14) 

r r 

Namely these conditions (B.14) fix an ambiguity in the sign choice in the relations (B.10) and 
(B.ll). Let us notice also that (B.l) in the limit y^ = y p = 0 takes the form 

Y _ N X 


sin (P±) = ±- 


cos(/3±) = ±- 


(B.15) 


Vl 2 + F 2 ’ V ' v/X 2 + T 2 ' 

By using (B.14) and (B.15) it is easy to see that two solutions ( a + ,(3 + ) and (a_,/3_) describe 
two opposite points of a unit sphere. As we explained earlier these two solutions describe two 
gyraton trajectories passing through the same space point (X, Y, Z). 

Using relations (B.5) and (B.6) one gets 

da(k = £,, d( k = |£| da = \/ r 2 - y 2 da . 

Let us denote 

i r 

(■ • -)y ± = 2 y o da sin a 5( dJk - C k) [■ • -]y p , 

then one has 

^ = -ZYp - A [ sin “ + <"' ,+ +*““-(•-ria^ • 


(B.16) 

(B.17) 

(B.18) 
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